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ABSTRACT 

By employing the equations of mean-square vorticity (enstrophy) fluctuations in strong shear flows, 
we demonstrate that unlike energy production of turbulent vorticity in non-rotating shear flows, the tur- 
bulent vorticity of weak convection in Keplerian disks cannot gain energy from vortex stretching/tilting 
by background shear unless the associated Reynolds stresses are negative. This is because the epicyclic 
motion is an energy sink of the radial component of mean-square turbulent vorticity in Keplerian disks 
when Reynolds stresses are positive. Consequently, weak convection cannot be self-sustained in Keple- 
rian flows. This agrees with the results implied from the equations of mean-square velocity fluctuations 
in strong shear flows. Our analysis also sheds light on the explanation of the simulation result in which 
positive kinetic helicity is produced by the Balbus-Hawley instability in a vertically stratified Keplerian 
disk. We also comment on the possibility of outward angular momentum transport by strong convection 
based on azimuthal pressure perturbations and directions of energy cascade. 

Subject headings: accretion, accretion disks — convection — hydrodynamics — MHD — turbulence 



1. INTRODUCTION 

Turbulence is usually driven by waves or instabilities in 
fluids, taping accessible free energy into fluctuating veloc- 
ity fields. For astronomical objects, accessible free energy 
can come from body forces {e.g., gravity or Lorentz force) 
or from large-scale fluid motions {e.g., rotation or shear). 

In Keplerian disks, thermal convection is an important 
ingredient in the thermal-viscous instability model which 
explain s the semi-periodic changes of ligh t curves of dwarf 
novae ( Meyer fc Meyer-Hofmeister 1981 ) and soft X-ray 
transients (Mineshige & Wheeler 1989). However, it has 
been argued that convective instability cannot access free 
thermal energy {i.e., higher entropy near the mid-plane) 
to transport energy vertically outward if it can access free 
rotational energy {i.e., higher angular velocities at small 
radii) to transport angular momentum radially outward. 
Recently a more general argument has been established 
to include the turbulence which is not thermally driven 
and whose azimuthal pressure perturbation is small. The 
equations of averaged turbulent kinetic energy due to lo- 
cal hydrodyn amic mixing can be expressed as (Balbus & 
Hawley 1998|) 
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where "losses" represent the energy sink due to viscosity. 



Vi is the turbulent velocity field, 5P is the pressure pertur- 
bation, p is the mass density, O is the Keplerian angular ve- 
locity, and the epicyclic frequency — {2QJr)d{r'^Vl) / dr . 
The notation < > in equations and (g) denotes the 
averaging carried out over the ensemble of turbulent cells, 
indicating the time/space correlation between two fiuctu- 
ating quantities. In the case of accretion disks, the time 
average is taken over several eddy turnover time (but still 
smaller than the viscous time scale), and the space aver- 
age is performed over 2tt in the azimuthal direction, over 
the whole disk vertical scale height^ , and over several eddy 
sizes in the radial direction (but still much smaller than 
the disk radius). 

As indicated by equation (^), the epicyclic term is an 
energy sink of {yl) if {vrVg) > 0. Therefore any local 
hydrodynamic instabilities cannot grow when azimuthal 
pressure perturbation is small. If a disk, however, is 
heated by other sources, the energy of weak convection 
can be maintained and then transport angular momentum 
inward. When thermal convection is weak, the associated 
negative Reynolds stress is a result of conservation of an- 
gular momentum of turbulent elements during the process 
of loc al mixing in R ayleigh-stable disks such as Keplerian 
flows ( [Balbus 2000| ). 

Besides considering the epicyclic effect which damps 
fluctuating velocity fields, the same damping effect should 
apply to fluctuating vorticity fields which characterize 
strong turbulence. It is also worth investigating how fiuc- 
tuating vorticity fields interact with background vorticity 
when convection is maintained by other heating sources 
such as M HD turbulence driven by the Balbus-Hawley 
instability (Balbus & Hawley 199S), in order to deter- 
mine th e typical linear modes for convec tion in nonlinear 
regime ( G!u, Vishniac, fc Cannizzo 200C| , hereafter GVC). 



lAirpraging rarnr tVip Hialf lipigVit is not performed whcn one would like to investigate the detailed vertical structures (e.g., see Stone e: 



al. 199e 



Miller fc Stone 2000.) 
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As high levels of fluctuating vorticities are observed in non- 
rotating shear flows at high Reynolds numbers, we must 
be wondering why the generation of fluctuating vorticities 
in Rayleigh-stable disks cannot be achieved by the usual 
vortex dynamics such as vortex stretchi ng or vortex tilting, 
as no ted and numerically simulated by Hawley, Balbus, & 
Wintc |rs 1999| . 

Perhaps one of the most difficult tasks in dealing with 
equations of turbulence is to determine the signs of cor- 
relations between fluctuating quantities. Positive, nega- 
tive, or no correlations are usually related to the properties 
of background flows such as stratification, rotation, back- 
ground shear, or body forces. Positive Reynolds stresses in 
non-rotating shear flows, for example, are always related 
to the mean flow moving in +x direction with negative 
gradient in +y direction. This situation is reversed in Ke- 



plerian flows for weak convection (Ryu & Goodman 1992; 


Kley et al. 1993; 


Kumcr et al. 1995;|Cabot 1996; Stone & 


Balbrjs 199^; pal 


3US 200C ) . Positive kinetic helicity occur- 



ring in magnetized disks shown in t he numerical simula- 
tion by Brandenburg & Donner 1997 probably means that 
the Coriolis force is of less importance in Keplerian disks 
where MHD turbulence is driven by the Balbus-Hawley in- 
stability. Similar to these examples, we would expect that 
background shear and rotation can determine the signs of 
correlations in the equations of fluctuating vorticities. 

Beginning with the perturbed vorticity equations in a 
non-rotating flow, we determine the signs of turbulent 
correlation based on positive Reynolds stresses without 
worrying about nonlinear shear instabilities in §2. In §3, 
we apply the analysis presented in §2 to Keplerian flows, 
showing that mean-square vorticity perturbations of weak 
convection cannot grow with positive Reynolds stresses. 
Positive kinetic helicity driven by the Balbus-Hawley in- 
stability in a stratifled Keplerian disk is analyzed through 
vortex equations in §4. In the last section, we comment on 
the possibility of outward angular momentum by strong 
convection in accretion disks based on the recent devel- 
opment of theories concerning about azimuthal pressure 
perturbations and directions of energy cascade. 

2. MOMENTUM CONSERVATION AND VORTICITY 
EQUATIONS IN NON-ROTATING SHEAR FLOWS 



The vorticity equation reads 



= (w • V) u— (V ■ u) w- 



Vp X VP 



Vxi/V^u, (3) 



where w = V x u, D/Dt = 9^-1- (u-V), P is the pressure, p 
is the mass density, and v is the kinematic viscosity. Con- 
sider a 3-D turbulent flow with a background shear dV/ dx 
without rotation, where V{x) is the background flow in +y 
direction. Let the vertical scale height be much smaller 
than the scale heights in x and y directions. If the eddy 
growth rate is smaller than the shearing rate {i.e., weak 
convection), linearizing the vorticity equation (u = V^J+v, 
and = V X v) gives the equations of mean-square vor- 
ticity (enstrophy) perturbation: 
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losses, (5) 
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where the terms "losses" represent the energy sink due to 
viscosity, Vi and uji are fluctuating velocity and vorticity 
respectively, and D/Dt denotes dt + Vdy. We have ignored 
the term (uiUjdjVi) which represents stochastic stretch- 
ing of vortex by turbulent shear. This effect is smaller 
than stretching by background shear as long as background 
shear rate is larger than turbulent growth rate. Similar to 
the role played by dSP/dr in equation (|l|), the terms as- 
sociated with dSP/dx and dSp/dx are energy sources of 
(i^y) in equation (^). Since these terms come from the 
baroclinic term in equation , they appear as a result of 
vertical stratification. 

The production or annihilation of turbulent vorticity re- 
lies on pressure perturbations, density perturbations, and 
the signs of correlations between fluctuating vorticities and 
fluctuating velocity shears as shown in the above equa- 
tions. However, the term associated with dy is usually 
small for weak convection as a result of strong shear. As- 
suming that all perturbation quantities are proportional to 
expiikxX + ikyy + ikzZ + iwt), we have the linear perturba- 
tion equations for adiabatic convection under Boussinesq 
approximation: 



ILOVy ■ 



dV 
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Vrr. 
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ituS — Vzdz In 



ik^^ 



iky^ 



pi/r 
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0, 
0, 
0, 

0, 
0, 



with the shearing constraint (Shu 1974) 



lUJ 



\dVldx\ 



(7) 
(8) 
(9) 

(10) 

(11) 
(12) 



in Shu 197. 



^This constraint is called the 
the shearing rate in this paper. In order to avoid confusion 



where Co is the frequency measured by a local observer co- 
moving with the mean flow, 5 = bpj p is the local fractional 
density perturbation, ^' = 5P/ p is the pressure perturba- 
tion divided by the density, Vx, Vy, and Vz are velocity 
perturbations, kx^ ky and kz are the wavenumbers in x, y 
and z directions. The validity of the plane-wave approx- 
imation for kx is assured by the shearing constraint for 

In fact, we are dealing with convection with the growth rates smaller than 
not use the term "weak-shear" here. 
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thermal convection in the presence of a shear flow. The 
plane-wave approximation for in the presence of a ver- 
tical gradient is not a bad assumption for weak convec- 
tion since first of all, the dispersion relation ( |l9| ) inferred 
from the above Boussinesq equations in the cas e of ac 



cretio n disks is the sam e as the one derived by Ruden, 



Papal aizou, fc Lin 198S , who consider the vertical wave- 
length of axisymmetrical modes for convection as a func- 
tion of z owing to vertical stratification. Second of all, as 
we shall see in the following text, the equations (^-(^2|) 
serve only as an example to illustrate how the phenomeno- 
logical interpretation works for explaini ng the directio ns of 
angular momentum transport done by Balbus 200C , and 
how convective eddies approximately look like in a shear 
fiow, all of which are not very sensitive to the detailed 
vertical structures of disks. We are going to adopt the 
same phenomenological strategy to evaluate the signs of 
correlations in the equations of perturbed enstrophy, and 
consequently to investigate if turbulent vorticities can be 
generated by vortex stretching/tilting in the regime where 
the terms associated with dy are small. 

By virtue of equations (Q) through (12) in the weak- 
convection limit (i.e., the convective growth rate iuj < 
\dV/dx\), one can show that the terms associated with 
dy (or say kg) in the above Boussinesq equations and in 
the equations of perturbed enstrophy are smaller than the 
terms associated with dV/dx, owing to the shearing con- 
straint. Without the pressure perturbation in the y di- 
rection, equation (|^) states that a turbulent element can 
conserve momentum in the y direction when it moves. 

In the case dV/dx < 0, conservation of momentum of 
turbulent elements in the y direction described by equa- 
tion (H) implies that {v^Vy) > 0, i.e., as observed in the 
comoving frame, positive (negative) values ofvx should oc- 
cur more frequently than negative (positive) ones when Vy 
is positive (negative), or vice versa. By virtue of equation 
®' iyl) ''^^ grow when {v^Vy) > ( Balbus fc Hawley 
1998|] The major concern for a successful growth of tur- 
bulent vorticity is that the first term of the right hand side 
in equation (|j) must be positive because the other terms 
associated with dy are small. Since it is the equation ^ 
which determines the sign of (vxVy) and the growth of tur- 
bulence, we would expect that the same equation plays the 
same role in equation (^). By virtue of equation (g|) again, 
we have 



2\dz dz 




where the growth rate 7 = i(I> > 0. The first term on the 
left hand is positive since turbulent velocity fields usually 
increase with convective growth rate. As a result, the cor- 
relation (dzVydzVx) is positive. The first term on the right 
hand side of equation (^ is « — (dzVydzVx) dV/dx which 
is therefore positive, leading to successful growth of (cj^). 
In fact, this positive sign can be imderstood in another 



way. I'i'he equation ot vorticity Uojx/Ut « clzVx(dV /dx) 
suggests that ojx tends to be negative (positive) when 
dzVx > (< 0) and dV/dx < (> 0). Hence lOx and dzVx 
are negatively correlated in a non-rotating shear flow. In 
other words, (w^) can grow from the interaction between 
background (dV/dx) and turbulent (dzVx) shear. 

Although the signs of correlations associated with back- 



ground shear in equation (|^) are not important because 
there are energy source terms due to pressure/density per- 
turbations, we can possibly determine the sign of (uJx'jJy) 
phenomenologically. In the case of dV/dx < 0, if we con- 
sider a turbulent element moving toward +x direction with 
positive sense of iOy, the element will tend to move toward 
+y direction more frequently so that the negative sense of 
LOx is created as seen in the comoving frame. This effect 
of vortex rotation is therefore described by the correlation 
(uJxi^y) < 0. The signs is reversed when dV/dx > 0. After 
all, the term (ojxOJy) dV/dx is an energy source of (ujy). 

The term (vxi-^z) in equation (^ describes the vorticity 
transport in the flows which have nonuniform background 
vorticity {i.e., d^V/dx^ ^ 0). Assume that the angular 
momentum of vertical vortex tubes are nearly conserved; 
i.e., the loss term in equation (p[) is small. If the gradient 
of background vorticity (P'Vj is positive (negative), the 
local turbulent mixing will transport the vorticity toward 
—X (-l-x) direction, or vice versa. This means that the term 
— {vxUJz) d'^V/dx'^ > 0. However, the production provided 
by vorticity transport is supposed to be small since Vx 
and the correlation length across the mean shear flows are 
severely shorten by strong shear. This point can be also 
understood in terms of the fact that the term (vxLOz) is 
zero in a shear fiow without a gradient of a background 
state in the x direction. Introducing a gradient of back- 
ground state gives rise to a factor about one over the scale 
height along the x direction which is much smaller than 
one over the vertical scale height associated with other 
terms in equation (|^). 

The second term in equation (^ describes vortex gener- 
ation due to the interaction between background vorticity 
dV/dx and variation of cross sections of vortices. The 
equation of motion DuJz/Dt k, —{dxVx + dyVy){dV/dx) 
suggests that uJz tends to be negative (positive) when 
dxVx + dyVy < (> 0) and dV/dx < (> 0). This 
means that the second term in equation (p) is the energy 
source of {<^z}- For a nearly incompressible disturbance, 
dxVx + dyVy < means vortex stretching {i.e., dzVz > 0), 
and dxVx + dyVy > means vortex squeezing {i.e., dzVz < 
0) in the vertical direction. 

After all, the turbulent vorticities lUx and LUy can grow 
simultaneously via vortex tilting by background shear 
{dV/dx) and by turbulent shear {djVi). The turbulent vor- 
ticity UJz can grow via vortex stretching/squeezing. The 
different vortex dynamics associated with uJz comes from 
the fact that lOz does not couple with lOx and ujy as shown 
in equations (^) , (|^) , and (^) . We note that vortex stretch- 
ing/tilting is a 3-D effect and cannot exist in 2-D fiuids 
{ujx ~ Wj, ~ 0) which are usually employed to model vor- 
tex generation in larger scales >> vertical scale height h 
{i.e., the shallow water approximation). Large-scale vor- 
tex production in 2-D usually relies on the baroclinic term 
in the first place, and subsequently vortexes can evolve by 
mutual interactions s uch as merging into large o nes which 
have long lifetimes (Adams & Watkins 199?;; podon & 

this evolu- 



Livio 1999). In this paper, we do not consider t 



tion. 

Before we apply the same approach to Keplerian flows, 
we need to emphasize that the above analysis just serves 
as preliminary calculations for Keplerian flows. The tur- 
bulent fields driven by thermal convection in a planar 
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shear flow are more complicated than those in a Keplerian 
flow. In the former case, thermal convection excites a non- 
linear sh ear instability which overwhelms the convective 
motion ( Stone fc Balbus 1996 ), and therefore the above 
analysis for thermal convection in a planar shear flow is 
quite questionable; for example, the te rm (ajjUJidiVi) which 



we ignored becomes most important ( Tcnnekcs fc Lumlcy| 
1972|)| In the case of Keplerian flows, however, numerical 



simulaitions thus tar have shown that strong epicyclic mo 



the azimuthal perturbed equation of motion which is sim- 
ilar to equation (||) except dV/dx is replaced by /2Vl\ 
i.e., 

"fvg = -—Vr = - Ir— +2n j Vr, (17) 

where we have dropped the term ikg^! . The above e qua- 
tion indicates that { vD cannot grow when (vrVg) > (Bal 



ipprcririctj iioii-ll ueai' riheat' liis tablllLieH (|Hawkiy, Bil^ 



tion s 

bus, ^ Winters 1999, however see Richard & Zahn 1999 



bus & Hawley 1998), which is different from the situation 



and Klahr 20001 ). Hence the above analysis should be rea- 



sonable for weak convection in a Keplerian flow. We shall 
see in the next section that the result implied from the 
equations of mean-square vorticity perturbations are in- 
deed consistent with the result inferred from the equations 
of mean-square velocity fluctuations. 

3. ANGULAR MOMENTUM CONSERVATION AND 
VORTICITY EQUATIONS IN A KEPLERIAN FLOW 

In Keplerian flows, the vorticity equations in cylindrical 
coordinates become 



Dt \ 2 / \ dz 
1 dp I dSP\ 
dz\ ^ rdO I 



K 

m 

losses, 



(14) 



D 
Dt 



dve 
dz 



^ J_9p / d6P 

dz \ dr 



K 

2h 



losses, (15) 



dvr 
dr 



dvg 
rde 



dr 



K 



2n 



losses, 



(16) 



where D/Dt denotes dt + ^dg. In a Keplerian disk 
threaded with sub-thermal magnetic fields, the "loss" 
terms in above equations include the turbulent damping 
due to radial mixing driven by the Balbus-Hawley instabil- 
ity. We do not show the terms associated with density per- 
turbations explicitly because unlike 6P/P <^ Sp/p in non- 
rotating shear fiows owing to equation (^, SP/P ^ Sp/p 
in rotating shear flows due to the Coriolis force (see equa- 
tion ( [l^ ) below). Once again ujz does not couple with ujr 
and LOg. 

Although perturbation equations (0) through (12) 



should be modifled by adding the Coriolis force in a rotat- 
ing flow, the Coriolis term 2il is comparable to the shear 
term rdfl/dr. Therefore the terms associated with dg are 
small compared to the first term on the right hand side in 
equation (p^, analogous to equation (Q). In a Keplerian 
flow, neglecting azimuthal pressure perturbations means 
that turbulent elements conserve angular momenta before 
they mix with background fluid. This leads to a transport 
of angular momentum down to the angular momentum 
gradient, givi ng rise to neg ative Reynolds stresses in a Ke- 
plerian disk (Balbus 2000). This picture is described by 



in a non-rotating shear flow owing to the appearance of 
the Coriolis force. 

In fact, negative Reynolds stresses can be realized as 
follows. Consider a turbulent element is moving with pos- 
itive Vr- While shear (rdfl/dr) tries to move the element 
to the +9, as observed in a corotating frame, the Corio- 
lis force (2ri) tried to move the element to the —9. Since 
20 > r\dfl/dr\ in Keplerian disks {i.e., > 0), the Cori- 
olis force is the winner. Therefore the turbulent element 
tends to have —vg, leading to {v^vg) < 0. Without rota- 
tion (Coriolis effect), (vxVy) > in a non-rotating shear 
flow with negative dV/dx. The contribution to {vrVg) from 
the radial deviation of turbulent azimuthal motion due to 
the Coriolis force is small since for weak convection in a 
Keplerian disk, equation (0) becomes 



ikr^ ^ 2Vlvg, 



(18) 



where the term iujVr is small and has been ignored. This 
means that the strong radial gradient of pressure pertur- 
bation is roughly balanced by the Coriolis force, resulting 
in negligible contribution to (vrVg). In other words, (vrVg) 
is mostly determined by equation (p7|). We note that in 
contrast to equation (ffl) , the extra term due to the Corio- 
lis force in equation (|^ indicates that the radial gradient 
of pressure perturbation should be large enough to make 
convection grow in a rotating flow. Namely, the Corio- 
lis force is a stabilized factor to convection. Owing to 
equations dl^) and (|l8|), the linear theory without radia- 
tive and turbulent damping indicates that 7^ is smaller 
than the Bnmt-Vaisala frequency —N'^ (GVC; Ruden et 
al. 1988): 



7 



A + l 



(19) 



where A = (kz/krY'. It has been known in terrestrial 
experiments that the critical Rayleigh numbers of rotat- 
ing Rayleigh-Benard conv ection are increase d above their 
non-rotating values {e.g., Zhong et al. 1993). 

The above analysis means that for weak convection, the 
signs of correlations in Keplerian flows in equations (14) 
(p^), and (|l^) are reversed compared with those in the 
non-rotating flows with dV/ dx < when we replace x — y 
with r — 9. In other words, we have (vrVg) < 0, {uirOJg) > 0, 
{uOrdzVr) > 0, and {ujz{drVr + dgvgjr)) < when k > 0. 
The right hand side of equation (|lj) becomes positive and 
{i^r) '^^^ grow. Although the first term on the right hand 



side of equation (15) turns out to be negative, (wg) can 



be maintained by the radial gradient of perturbed pressure 
owing to strong epicyclic effect. 

If we assume that weak convection transports angular 
momentum down to the angular velocity gradient, then the 
signs of correlations in equation (nih and equation (|5h are 
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the same as those in the non- rotating flows with dV/ dx < 
when we replace x — y with r — 9. This means that the first 
term on the right hand of equation ([l^ ) is negative because 
> 0. Hence (w^) cannot grow and weak convection dies 
away. 

If the disks are Rayleigh-unstable {i.e., < OV the 
signs of aU terms in equations (p^), and ( |l6|) re- 

main unchanged compared to those with positive ex- 
cept that the first term on the right hand side in equation 
( p^ ) becomes positive because {uJrUJe) < 0. In other words, 
a Rayleigh-unstable disk behaves as a non-rotating shear 
flow in the sense that (vrVe) > 0. We note that our analy- 
sis for Rayleigh-unstable disks is oversimplified because a 
non-linear shear instability should be excited. 

4. POSITIVE KINETIC HELICITY AS A RESULT OF THE 
BALBUS-HAWLEY INSTABILITY 

A weakly magnetized Keplerian disk is linearly unsta- 
ble to the Balbus-Hawley instability. In nonlinear regime, 
this instability can drive MHD turbulence which in turn 
leads to a dyna mo process unless the rn a gnetic Reynolds 
numbers are l ow ( Balbus fc Hawley 1998 ). Brandenburg & 
Donncr 1997 shows that this dynamo process observed in 
simulations can be imitated by an a — 17 dynamo in a verti- 
cally stratified disk, with a negative agg in the upper disk 
plane (where agg is the 9-9 component of kinetic helicity 
tensor a). A negative agg means that the kinetic helic- 
ity (vzLOz) is positive for nearly isotropic turbulence. This 
contradicts with the usual notion that the kinetic helicity 
of con vection driven by th e Coriolis force is usually neg- 
ative. Brandenburg 1998 suggests that the right-handed 
helical turbulence results from the combined effect of the 
Balbus-Hawley instabihty {i.e., (SBrSBg) < 0) and mag- 
netic buoyancy. 

Although we have concentrated on convection-like ed- 
dies and have not taken into account the effect of sub- 
thermal magnetic fields in this paper, equations (^) and 
( p^ ) could shed light on preferred directions of helical tur- 
bulence if we realize that a Coriolis-dominated disk char- 
acterized by equation (^6|) becomes a shear- dominated 
disk symbolized by equation (||) due to the Balbus- 
Hawley instability. As explained in the preceding sec- 
tion, {uJz{drVr + dgvg/r)) < in a nonmagnetized Kep- 
lerian disk. Assume that a turbulent element expands 
(contracts) as it goes upward (downward) in the upper 
disk plane. The above negative correlation means that a 
turbulent element tends to have a negative (positive) lUz 
due to expansion (contraction). Consequently, (vzCOz) < 
(see the vortex plot associated with positive in Figure 
1). This result is consistent with the picture in which the 
Coriolis force is the winner over shear so that turbulent 
eddies are mostly left-handed; i.e., as observed in a coro- 
tating frame, left-handed helical turbulent motion results 
from a positive vorticity (k^ > 0) of the mean flow. How- 
ever, introducing the Balbus-Hawley instability changes 
the sign of Reynolds stresses from negative to positive in 
Keplerian disks because angular momentum of a turbu- 
l ent element i s changed by fluctuating magnetic torques 
(Balbus 2000), leading to the condition for the instabil- 



ity changed from the Rayleigh criterion k < to the 



Chandrasekhar-Vclikov crit e rion dQ /dr < (Velikhov 



1959 



criterion is a manifestation of conservation of angular ve- 
locity of turbulent elements due to the fact that the fluctu- 
ating magnetic torque counter acts the torq ue of the Cori- 
olis force in a weak field limit ( Fricke 1969| ). Without the 



torque by the Coriolis force, azimuthal dynamics of tur- 
bulence in Keplerian disks becomes that in a non-rotating 
shear flow. Furthermore, if the contribution to {vrVg) due 
to the radial deviation of azimuthal motion is small (like 
the situation described by equation (18); radial deviation 
of azimuthal motion is usually small owing to shear), then 
the signs of correlation quantities in a weakly magnetized 
disk should behave as the one described by negative vor- 
ticity of the mean flow (such as the correlation quantities 
in equations (0), (^, and (^), giving rise to a positive ki- 
netic helicity (see the vortex picture associated with neg- 
ative dV/dx in Figure 1). The only physical difference of 
vortex dynamics between a non-rotating shear flow and a 
weakly magnetized disk is that in the later case, a pos- 
itive kinetic helicity coupled with a positive vorticity of 
the mean flow becomes an energy sink of {(^z) shown in 
equation ([l6|). This is analogous to the result that a pos- 
itive Reynolds stress is an energy sink of Vg in Keplerian 
disks ([Hawley, Balbus, fc Winters 1999)). 




y 




dV/dx > 



dV/dx < 



phandrasekhar 1960| ; Fricke 196S) in a Keplerian 



flow. Equivalently speaking, the Velikov-Chandrasekhar 



Fig. 1. The expansion of a rising turbulent element described 
by equation ( p^ or by equation (g) when > or dV/dx > 
(left), and when < or dV/dx < (right). As observed in 
a corotating frame, {LOz{drVr + dgve /r)} < and {v^uj^) < 
when K'^ > 0. These signs are reversed when < 0. A nega- 
tive (positive) ojz, as a result of a positive (negative) vorticity 
of the mean flow, can give rise to a negative (positive) (vrVe) 
when radial deviation of azimuthal motion is small. As argued 
in the text, the picture on the right could also apply to the 
uJz produced by the Balbus-Hawley instability which generates 
positive Reynolds stresses in Keplerian disks. In the magnetic 
case, however, the vortex motion illustrated in the figure be- 
comes an energy sink oi {u)1^. ke < kr is a, result of shear. 



Although our analysis is related to convection-like eddies 
{i.e., kz ^ l//i), the picture we present here for estimating 
the sign of kinetic helicity does not contradict the result 
that the Parker instability, a magnetohydrodynamic in- 
stability with kz ^ 1/h, is su ppressed by MHD turbulence 
( Vishniac fc Diamond 19*9^ ). The main point presented 
here depends solely on the sign of Reynolds stresses which 
is assumed to be largely caused by the azimuthal deviation 
of radial motion. Therefore, the picture that we sketch 
here applies to any turbulent eddies which expand (con- 
tract) as they rise (descent). We note that rising turbu- 
lent eddies can expand statistically even without any types 
of buoyancy. As a matter of fact, the vertical structures 
of large-scale magnetic flelds generated by the Balbus- 
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Hawley instability in vertically stratifi ed accretion disks 
are stable against m agnetic buoyancy ( ^tone et al. 1996 ; 
Miller & Stone 2000). The detailed information, such as 



the magnitude of this kinetic helicity and the expansion 
of turbulent eddies in a stably stratified disk, should de- 
pend strongly on the detailed structure of the flow, such 
as azimuthal perturbations {i.e., anisotropic turbulence) 
and correct vertical dynamics. 

5. DISCUSSION AND OUTLOOK: COMMENTS ON 
ANGULAR MOMENTUM TRANSPORT BY STRONG 
CONVECTION 

In t his paper, we demonstrate that weak convection can- 
not bi p self-sustained by vortex stretching and tilting in a 



and by MHD turbulent damping. From equations ([141), 
(|l|),and([l6l), we have seen how convective vorticity /shear 
possibly interacts with background vorticity/shear and 
how this interaction becomes an energy sink of turbulent 
vorticity. This means that the quasi-linear approach for 
the mixing-length theory by GVC based upon the sec- 
ondary Kelvin-Helmholtz instability is qualitatively cor- 
rect. 

However, the issue of angular momentum transport 
by weak convection is probably not of great importance 
in reality. According to the thermal-viscous instability 
model, weak convection occurs in K eplerian disks when 
the S hakura-Sunyaev viscosity ass ( Bhakura fc Sunyaev 
1973) is larg e ~ 0.1, a s a result of strong radiative 



Keplerian flow, suggested by the equations of enstrophy 
fluctuations. This is because epicyclic term is an energy 
sink of {uir} when {vrVg) > 0. The situation is reversed 
when weak convection has a negative Reynolds stress, lead- 
ing to a successful growth of thermal convection. Our de- 
termination of the signs of correlation quantities in the 
equations of enstrophy is based on the point of view in 
which the signs of Reynolds stresses are determined by 
conservation of angular momentum for weak convection 
( Balbus 2000 ) . Although our determination of signs of cor- 



losses (GVC; Gabot 1996 ). Furthermore, based on GVC, 
large ass means strong MHD turbulent mixing driven 
by the Balbus-Hawley instability, which suppresses weak 
convection by smoothing away momentum and entropy 
anisotropics associated with convective bubbles mostly in 
the radial direction (see equation (|9|) below): ^ << 1 
when — iV^ << fi). Consequently, only the convection 
with larger Rossby numbers Ro ^ 1 can survive in a mag- 
netized Keplerian disk. Weak convection can occur at 
smaller r adii in protostellar disks (e 



see 



relati( |)n quantities is not a rigorous proof but quite sugges- 
tive, the results are consistent with those suggested by the 
equations of velocity fluctuations. Full three-dimensional 
numerical simulations are needed to verify our thought. 

Positive kinetic helicity in a weakly magnetized, verti- 
cally stratified Keplerian disk impl ied from our analysis 



D'Alessio ei 



al. 1998). In such a cold disk, a significant 'dead zone 



agrees wi th the simulation result by Brandenburg & Don 



ner 1997. However, the physical picture shown in Figure 
1 is totally different from the plot sketched in Figure 7 
in their paper: in a corotating frame of a Keplerian disk, 
the direction of Coriolis effect (20.) should be opposite to 
the direction of shear effect (dil/dr), resulting in a pos- 
itive vorticity (k^ > 0). In our analysis for turbulence, 
{wz{drVr + deve/r)) < (> 0) when vorticity is posi- 
tive ( negative), giving rise to negative fpositive) (v-^uj-,) as 



implies that the Balbus-Hawley instability may not func- 
tion properly, but accretion might occur near the surfaces 
of disks where cosmic rays ( Gammie 1996| ), irradiation 
(Glassgold et al. 1997), or rigorous heating beyond the 
bottom of photosphere as manifested by protostellar winds 
and flares mig ht increase the ionization le vel beyond the 
critical values (Gammie and Menou 1998). However, the 
heat due to layer accretion should reduce the tempera- 
ture gradient of the 'dead zone' in the vertical direction, 
possibly leading to the suppression of weak convection as 
a result of enhanced radiative losses in the radial direc- 
tion (GVC). Therefore, the question becomes, does strong 
convection transport angular momentum i nward or out- 



ward? Lo w Reynolds number simulation by Cabot & Pol 



long 4s {vAdrVr + dovo/r)) is positive. Since the Balbus- lack 1992| found strong convection transports angular mo 



Hawley instability generates a positive (vrVg) which resem- 
bles the hydrodynamic turbulence with a negative vortic- 
ity, MHD turbulence driven by the Balbus-Hawley insta- 
bility should display right-handed helical motion as ob- 
served in numerical simulations. However, if rising eddies 
contract statistically, our analysis will give rise to the con- 
ventional sign of kinetic helicities in a rotating fluid in- 
fluenced mostly by the Coriolis force. The expansion of 
rising turbulent elements in a stably stratified disk could 
be just a result of a special profile of underlying back- 
ground states, or might result from the possible scenario 
that turbulent mixing, an effect definitely not included in 
usual laminar analyses, alters the condition of flux freezing 
in some degree. On t he contrary, the simple explanation 
by Brandenburg 199? relying on magnetic buoyancy is not 
consistent with the simulations for vertically stratified ac- 
cretion disks stirred by the Balbus-Hawley instability. 

GVC studied convective energy transport in Keplerian 
disks stirred by the Balbus-Hawley instability. The typical 
modes of convection are determined from linear perturba- 
tion theory and nonlinear saturation which is caused by 
the secondary Kelvin-Helmholtz instability (i.e., interac- 
tion between convective shear and background vorticity) 



mentum outward. 

Azimuthal pressure perturbations have drawn large at- 
tention in the literature recently, as suggested by the 
equations of mean-square velocity perturbation and by, 
in this paper, the equations of enstrophy. Axisymme- 
try/nonaxisymmetry of convection is actually an intimate 
issue of the damping of convection due to MHD turbulent 
mixing (|Klahr, Hcnning, fc Kley 1999|). When MHD tur- 



bulent mixing is large (i.e., large ass), modes with large 
kr are all suppressed owing to strong radial mixing. Sim- 
ilarly, modes with large kr should also succumb to strong 
radiative losses in radial direction. Modes with small kr 
can survive in a shearing environment only if kg is also 
very small, suggesting nearly axisymmetrical convection. 
As mentioned in the preceding paragraph, weak convection 
occurs when ass is large. Nearly axisymmetrical convec- 
tion is therefore a manifestation of weak convection which 
struggles with turbulent and radiative damping in a differ- 
entially rotating flow. In other words, nonaxisymmetry of 
strong {Ro ~ 1) convection means that azimuthal pressure 
perturbations could be important to violation of conserva- 
tion of angular momentum. If this happens, the equations 
of velocity and vorticity fluctuations could allow strong 
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convection to transport angular momentum radially out- 
ward. Without a doubt, the quasi-linear analysis based 
on equations of velocity and vorticity fluctuations is very 
suggestive. 

The other point of view regarding the nature of posi- 
tive/negative viscosity is the direction of turbulent energy 
cascade. Negative eddy viscosity manifests the process in 
which turbulent fields do not only extract energy from the 
mean flow, but the extracted energy is also passed up to 
larg er scales and finally goes back to the mean flow ( Starj 



necessarily true. Klahr 200C shows that circumstellar disks 
with outward transport of angular momentum by convec- 
tion cool continuously throughout the simulation. 

Direction of angular momentum transport by convec- 
tion is also an important issue in advection-dominated ac- 
cretion flows (ADAFs). When ass (not due to convec- 
tion) is small, strong convection [Ro 1) occurs when 
convection transports angular momentum inward (or out- 
ward less efficiently than it tra nsports energy), and ac- 



1968) The connection between negative viscosity and in- 



verse energy cascade in two-dimensional flows has been 



cretion is suppressed i n ADAFs (Narayan, fgumenshchev 
fc Abramowicz 2000| ; Quataert fc Gruzinov 2000| ). How- 



studied in fluid socie ty (Kraichnan 1976; Pouquct 1978 
Chechkin et al. 1999). Cabot 1996 and Klahr et al. sug- 



gested that negative viscosity is a result of an inverse en- 
ergy cascade of nearly axisymmetrical convection and that 
positive viscosity results from nonaxisymmetrical convec- 
tion in which 3D hydrodynamical turbulence is resumed. 
In fact, the collapse of 3D hydrodynamical turbulence to 
2D is not necessarily related only to axisymmetrical pat- 
terns, but is a natural result of turbulent anisotropy due to 
strong body forces, which has been manifested by ample 
examples in geophysics. Hossain 1994 found that when a 
strong rigid rotation [Ro < 1) is turned on, the turbulent 
velocity fields perpendicular to the direction of rotation 
{i.e. X — y plane) are strongly correlated along the direc- 
tion of rotation (i.e. 2; direction) . Therefore 3D turbulence 
(direct cascade) reduces to an approximate 2D state (in- 
verse cascade). In the case of accretion disks, epicyclic 
effect should enforce turbulent fields of weak convection 
to have an approximate 2D state, giving rise to negative 
viscosity due to an inverse energy cascade. If convection is 
strong, however, the collapse to 2D breaks down and this 
probably results in positive/zero eddy viscosity due to a 
direct energy cascade. 



ever, as noted in this section, strong convection is unlikely 
to be nearly axisymmetrical and is doubtful to be able 
to collapse to a 2D state especially in a thick disk where 
the direction of strong entropy gradient has a component 
along angular momentum gradient. Numerical observa- 
tions of a transition of direction of angular momentum (if 
it exists), a transition of direction of energy cascade, and 
a transition of pattern formation for convection in both 
Keplerian accretion disks and ADAFs is worthwhile when 
Ro and Ra are the control parameters, where Ra should 
be the "turbulent" Rayleigh number characterized by ra- 
diative/ ad vection losses and MHD turbulent damping. 

While the perturbed baroclinic terms in the equations 
of perturbed enstrophy act as the energy source of fluctu- 
ating vorticities, we do not consider the terms associated 
with the background baroclinic term in this paper, due to 
the thin-disk approximation. That is, Q does not vary 
with z dramatically in thin d isks, resulting in a barotropi c 
and Keplerian disk {e.g., see Frank, King, & Raine 1992). 



The numerical simulation by Klahr 2000 found that the 



gradient of fluctuating angular momentum generated by 
weak convection could flatten background gradient of an- 
gular momentum. Consequently, a secondary instability 
occurs at a fast rate of some fraction of Q. This induced 
instability overwhelms weak convection and drives angular 
momentum outward. In terms of the quasi-linear analy- 
sis based on velocity and vorticity fluctuations, outward 
angular momentum transport by this secondary instabil- 
ity is a result of strong azimuthal pressure perturbations 
and stochastic stretching of vortex by turbulent shear. By 
virtue of nonlinear phenomena of turbulent cascade, this 
fast secondary instability diminishes epicyclic effect and 
thereby forces fluid back to an approximate 3-D hydrody- 
namical state. 



If t lermal convection generates positive eddy viscosity 



The background baroclinic effect is usually of less impor- 
tance on local faster turbulence in disks; for instance, the 
Biermann battery owing to the same mechanism is usu- 
ally ignored when a mean-fleld dynamo is able to operate 
via mean kinetic helicities driven by local turbulent events 
such as thermal convection in stars or supernovae/stellar 
winds in galactic disks. In spite of their appearance in 
thick disks which behave as fast rotating stars, it is also 
quite doubtful that the background baroclinic term can 
have important influence on the onset of convective insta- 
bility in these cases. As long as the convective growth rate 
is not too small, the effect of large scale circulation due to 
the baroclinic term to smooth away entropy inhomogeneity 
associated with small-scale mixing such as thermal convec- 
tion in this case is usually negligible, compared with the 
same damping mechanism achieved by MHD turbulence. 
On the other hand, as we mention in the end of §2, a 
large-scale vortex can be initiated from the background 
baroclinic effect. Large-scale ed dies could be also formed 
via merging c onvective vortexes (Zhong et al. 1993; Godon 



in Keplerian disks, it does not necessarily mean that it 
can be self-sustained. It has been a concern that a self- 
sustained convection would violate the second law of ther- 
modynamics: in the case of a system in which the eddies 
are thermally driven, the heat dissipated fro m the mea n 
flow cannot again be used to drive the eddies ( Starr 1968| ) . 
If thermal convection can be self-maintained in Keplerian 
disks, this could mean that convection is not totally ther- 
mally driven, in the sense that it can access the available 
rotational energy a t smaller radii via s econdary or non- 
linear instabilities (Kumer et al. 1995; Richard fc Zahn| 
19990 However, the reverse of the above statement is not 



fc Livio 1999) or throu gh hclicity fluctuations of convection 
^Branover et al. 1999). The nonlinear interaction between 



large-scale vortexes and convective vorticities is worth in- 
vestigating. 

Finally we note that our analysis in this paper is applied 
to local turbulent mixing in accretion disks, such as ther- 
mal convection and the Balbus-Hawley instability. Turbu- 
lent transport mediated by global waves, such as internal, 
Rossby, or spiral shock waves, is not subject to our analy- 
sis. Moreover, convective turbulence might behave differ- 
ently in the environment of MHD turbulence since ideal 
invariants during turbulent cascades are different between 
hydrodynamical and MHD turbulence {e.g., sec Biskamp 
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1997). 
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